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Abstract 

We investigate the fast-reaction asymptotics for a one-dimensional reaction- 
diffusion (RD) system describing the penetration of the carbonation reaction in 
concrete. The technique of matched-asymptotics is used to show that the RD 
system leads to two distinct classes of sharp-interface models, that correspond 
to different scalings in a small parameter e representing the fast-reaction. Here e 
is the ratio between the characteristic scale of the diffusion of the fastest species 
and the one of the carbonation reaction. We explore three conceptually different 
scaling regimes (in terms of e) of the effective diffusivities of the driving chemical 
species. The limiting models include one-phase and two-phase generalised Stefan 
moving-boundary problems as well as a nonstandard two-scale (micro-macro) 
moving-boundary problem - the main result of the paper. Numerical results, 
supporting the asymptotics, illustrate the behavior of the concentration profiles 
for relevant parameter regimes. 

Keywords: Concrete carbonation; Reaction layer analysis; Matched asymptotics; 
Fast-reaction asymptotics; two-scale sharp-interface models; numerical approxima- 
tion of reaction fronts 

Abbreviations: FBP - free boundary problem; PDE - partial differential equation 

1 Introduction 

Carbonation reactions alter the concrete's pore geometry in a difficuh-to-control 
fashion. A good theoretical multiscale understanding of the evolution of carbonation 
reactions in such a complex multiphase material is vital towards obtaining accurate 
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predictions of the durability of large concrete structures (see, for example, \^). 
Having a good estimate on the life service of motorways, bridges, sewage pipe systems 
(etc.) can save significant amounts of money and energy yearly - hence the growing 
multidisciplinary research interest in this topic. 

1.1 Background 

From the point of view of the involved geochemistry, the process of concrete car- 
bonation mainly involves the reaction of atmospheric carbon dioxide with calcium 
hydroxide found in the pore solution to form water and calcium carbonate. The pore 
solution is highly alkaline (pH ~ 14) and as soon as the pH level decreases the micro- 
scopic oxide layer at the steel reinforcements disappears and the steel can corrode. 
This corrosion causes the durability of the structure to reduce dramatically. Figure 
[T] shows how the pH levels of the concrete block drop dramatically. For more details 
on concrete carbonation and its relevance with respect to corrosion and durability 
issues, we refer the reader to [2H], [IS] and |32j as well as to the references cited 
therein. 




Figure 1: Cross-section of a reinforced concrete block. Pink indicates low pH while 
colourless indicates a high one. 

Here, we focus on the dominant carbonation reaction mechanism, namely 

C02{g aq) + Ca(0H)2(s ^ aq) Ca.CO^{aq ^ s) + H2O. (1) 

In other words, atmospheric carbon dioxide diffuses through the unsaturated con- 
crete matrix and dissolves into the pore water while the calcium hydroxide is avail- 
able in the pore solution by dissolution from the solid matrix. Free water and calcium 
carbonate are the main reaction products in ([T|). CaCOs quickly precipitates into 
the solid matrix, while the water diffuses through the material very much influenced 
by eventual wetting and deleting cycles. As a consequence of so many reaction and 
transport mechanisms, the available macroscopic models capturing carbonation in 
concrete are complex and involve a large number of parameters. For this reason, al- 
though the mechanisms are rather well understood, the predictability power of such 
models is questionable. The challenge is to find minimal models that can accurately 
predict CO2 penetration depths in concrete. 
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Our objective is to understand how fast the reaction ([T]) moves into the material, 
driven by CO2 diffusion. In other words, we want to use a detailed reaction layer 
analysis to determine the speed of growth of the macroscopic carbonated phase (the 
pink region in Figure [T]). The output of our approach/models consists of carbonation 
penetration depths and profiles of active concentrations. 
In this context, the open question is: 

(Q) Which boundary conditions need to be imposed at the reaction interface to 

describe its correct motion? 

(Q) has been recently addressed, for instance, in [2(H IMj [2| [23| [22| [3| [2T] . and [25] . 

where the authors indicate that the answer to (Q) seem to be rather well under- 
stood in one-space dimension. However, the 2D and 3-D cases are comparatively 
untouched. The reason is simple: For the moment, one does not know how the re- 
action front [hosting ([1])] feels the presence of corners or more complex geometrical 
shapes. Furthermore, the mathematical theory of FBPs in more space dimensions 
[and posed in heterogeneous media] is less understood. Rigorous results are mostly 
available for standard FBPs, relying heavily on a lot of a priori prescribed regularity 
for the moving boundary (sharp interface). For practical purposes, we need robust 
models [e-approximations of sharp-interface models for carbonation] that are appro- 
priate in any space dimension. Such models should not require any special regularity 
requirements on the sharp interface and should be able to capture as e — )• the cor- 
rect moving interface conditions. This is the place where our paper contributes. The 
research reported here is preliminary, with our main asymptotic results emphasizing 
a new, non-standard two-scale sharp-interface model. Many associated fundamental 
and practical issues are currently open (see section [8]) and deserve further study. 

1.2 Approach 

We introduce a basic macroscopic reaction-diffusion model to describe the aqueous 
chemistry and transport involved in ([T]). To be more precise, we only consider the 
concentration of the reactants C02(aq) and Ca(0II)2(aq) which we allow both to 
diffuse. Furthermore, we assume that the diffusion coefficients can be concentration 
dependent and also may vary with the reaction rate. We then use the technique 
of matched-asymptotics [5] to perform a detailed analysis of the reaction layer for 
various physically relevant scaling regimes. Additionally, we use COMSOL to pro- 
duce numerical results of the original system that corroborate the asymptotic results 
obtained for the various scaling cases. 

Since we look to carbonation in high Thiele moduli regimes, the rate at which the 
concentrations diffuse is typically much slower than that of the reaction. Thus we 
focus on the fast reaction limit e — t- 0, where e > is a scale parameter that will be 
explained in the next section. Assuming strong time scales separation, the reaction 
is expected to be in equilibrium except for an interior narrow layeiQ where virtually 
all the reaction takes place. The structure of this localized reaction zone is sensitive 

^Interestingly, the strong separation of the characteristic time scales of diffusion and reaction 
lead to the presence of two distinct (macroscopic) characteristic space scales that relate to the 
typical sizes of the layers hosting diffusion and reaction. 
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to the behaviour of the diffusion coefficients. For this reason, we investigate three 
different possible situations in regard to the diffusivities. These lead to the derivation 
of conceptually different sharp-interface models. We see three main types; namely, 
one and two phase Stefan like problems with zero latent heat and two scale micro- 
macro free boundary problems. The latter is a completely new type of model in 
which the speed of the reaction interface is updated from a smaller length scale. 

The broad spectrum of different models showed here depicts not only the be- 
haviour of the concentrations for many possible scenarios of carbonation, but they 
can also be related to other types of aggressive reactions like redox scenarios |14j . 
silicon-oxidation , sulfate attack [H [121 ISS] and combustion [13] . The techniques 
used here are therefore applicable to a larger set of reaction-diffusion problems. 

1.3 Organization of the paper 

The macroscopic reaction-diffusion system modelling carbonation is introduced in 
Section [2j Performing our fast-reaction asymptotics, a first sharp-interface model 
is derived in Section [3l The bulk of the paper (Sections [H [5] and [6|) is devoted to 
the reaction layer analysis for three different diffusion coefficient regimes. The first 
regime (Section S]) is that in which the diffusion coefficients remain order 1, which 
we refer to as slowly varying. This is the usual case which has received attention 
in the literature. The other regimes, allow the diffusion coefficients to change their 
orders of magnitude in terms of the parameter e, the diffusion coefficients then 
being termed rapidly varying. Thus, as a second regime (Section [5|) we consider 
only a rapidly varying CO2 diffusivity, with the hydroxide diffusivity slowly varying. 
Finally, as a third regime (Section [6]), we consider the effects of rapidly varying CO2 
and hydroxide diffusivities simultaneously. The main results of the paper are the two- 
scale sharp interface models, which are summarised in [71 The specific generalised 
Stefan and kinetic conditions on the moving boundary are summarised for the single 
and two rapidly varying cases in sections [5.31 and [U^ and should be compared to the 
conditions in the slowly varying case in section [^31 Each regime subdivides according 
to the size of a relative transport parameter, which measures the characteristic ratio 
of diffusivities and concentrations for CO2 and hydroxide. 

We close the paper with a discussion section around possible connections between 
the nonstandard two-scale models derived in this paper with the two-scale models 
obtained previously in the homogenization literature |15j . 

2 The model equations 

We consider the concrete occupying a one-dimensional slab geometry :=]0,L[. 
The X-axis be directed into the concrete with the surface x = being exposed to 
an external source of CO2, whilst the surface x = L is assumed impervious to all 
reaction species. We let c = c{x, t) and h = h{x, t) denote the concentrations of 
C02(^qj and Ca(0H)2(aq) respectively, expressed as moles per unit volume (i.e. the 
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intrinsic concentrations). We adopt the following set of reaction-diffusion equations, 
InO<x<L,t>0 



with boundary conditions 

onx = 
onx = L 
and initial conditions 

att = 

The initial hydroxide concentration has a representative value (taken as maximum) 
denoted by the constant h^, the most common situation being /ij(x) = l,Ci(x) = 
for < a; < L. Here the concrete substrate has length L; R is the carbonation rate 
with reaction rate coefficient k; Dc and Dh are the diffusion coefficients of CO2 and 
Ca(0H)2 respectively, which are not necessarily constant. Not only can the diffusion 
coefficients be concentration dependenlH but we also allow them to vary with the 
reaction rate. Motivation for this comes from the significant change that takes in 
the concrete matrix during carbonation. The main modelling assumption is thus to 
consider their dependence linked to R{c,h), which for specificity may be taken in 
Arrhenius forms 

Vc, T^h suitable constants and Rc, Rh are based on partial reaction rates. We note 
that neither the effective diffusivities nor the reaction rate coefficient depend on 
humidity. Also the transfer of CO2 from the air to water phase (and vice versa) and 
the dissolution of Ca(0H)2 from the solid matrix to water phase (and viceversa) 
are in local equilibrium i.e. all production terms by Henry's law vanish; see [6] for 
more details on molecular transfer across water-air interfaces and on Henry's law. As 
such, from a modelling perspective, we may justify the Robin boundary condition 
for the aqueous carbon dioxide in Q as follows. Denoting the concentration of 
gaseous carbon dioxide C02i^^^y in the concrete by Cg with diffusion coefficient Dc^-, 
its exchange with the external atmospheric concentration c* at the exposed concrete 
surface is given by 

^Consider the diffusion coefficient of CO2, this may be expressed as Dc = -DcW$(/i), where w is 
a tortuosity factor and is the porosity which depends on the concentration of Ca(Oh)2. This 

gives a justification for the CaOH2 dependance of the diffusivity. We refer the reader to [5D] for 
more remarks in this direction. 



dc d / dc 

dt dx \ ^ dx 

dh _ d_f dh 

dt dx \ dx 



R{c,h), 
-R{c,h)- 



(2) 
(3) 



- DcTT = H [c - c), — = 0; 
dx dx 

l^ = 0'? = 0^ 
ox ox 



(4) 
(5) 



c = c*Ci{x), h = hPhi{x). 



(6) 
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where H* is the mass transfer constant. 

Assuming an equihbrium balance c = C^Cg for the local exchange of gaseous and 
aqueous carbon dioxide, we may combine these expressions to give ([U, where 



TT* _ 9 ^ * _ r<H * 

H - C -C Cg 



The dimensionless Henry constant is temperature dependent, a typical value 
being 0.82 at 20°C. We remark that the assumption of equilibrium balance between 
gaseous and aqueous carbon dioxide means that Dc should be taken as Z?cgi the 
governing equation ([2]) then being consistent with the models of [29] when written 
in terms of Cg. More sophisticated models may relax this equilibrium assumption, 
allowing CO2 in the gas and aqueous phases to be considered separately. 

The formation of calcium carbonate (concentration z{x, t)) can be modelled using 
the rate equation 

dz 

g^ = i?(c,/i), (8) 

where its diffusivity is taken as being negligible. As such the amount of carbonate 
can be determined once ([2])^® is solved together with specifying a suitable initial 
condition e.g. z = at t = 0. 

We non-dimensionalize as follows, 

^2^0 _ 

X = Lx, t = — 7^ — t, c = c*c, h = h?h, 
D^c* 

Dc = D^cDc, Dh = DlDh, R{c, h) = 9r{c, h), 

using L as the characteristic length scale and 6 a representative reaction rate scaling. 
We have taken D^ and 
Dropping" 's, we obtain 



We have taken D^ and D^ as the maximum values of the diffusion coefficients. 



In0<x<l,i>0 e'U — Dc^]]=-ric,h), (9) 



with boundary conditions 



dt dx \ dx J 
•h _ 6^d_ / dh 
dt dx \ ^ dx J 



n 1' dh S"^ d / dh\\ , , , , , 



onx = -Dc^ = H{l-c), ^ = 0, (11) 

onx^l |^0,|^0; (12) 
and initial conditions 

att = c = Cj(x), h = hi{x), (13) 

where 

,2_D'cC* .2_Dlh^ ,^_c* LH* 

'^"W ^-w ^^^^ 
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In the semi-infinite concrete case, L is at our disposal and can be taken to be D^/ H* 
so that H = 1. 

The scaling 9 is chosen so that r(c, h) < 1 for 0<c<l,0</i<l. Tobe more 
precise, for reaction terms of the form 

R{c, h) = kd'h'i 

where A: is a constant reaction coefficient and p and q are positive constants, we have 

r{c,h) = d'K' (15) 

by choosing 6 = k{c*Y{hPY. Referring to [31] and [19], typical values of the di- 
mensional parameters are given in Table [TJ A representative length L can be of the 
order of cms upto several metres i.e. L = 0.1 — 10m depending upon the situation 
and geometry considered. The reaction rate corresponds to the common partial re- 
action orders p = q = 1. However, other values of these, 0.9 < p,q < 1.5, also seem 
appropriate. These give the estimates showned in Table [2j 
The parameter range of relevance is thus 

e<l, <5<0(e), fi^l, H = 0{1), (16) 

where we note that the "relative transport parameter" 

- = (17) 

is typically small or at most order 1 [28]. The CO2 interfacial exchange parameter H 
is estimated using the maximum values of the CO2 diffusivities, which are expected 
to occur at the surface and for gas are denoted by . Since Z)f? = D^. we have that 

O J Cg C Cg 

H = LH* / D^^ , the larger values of this parameter suggesting that the interfacial 
exchanges are close to equilibrium. The diffusion coefficients now have dependencies 
on the parameters e and 5 i.e. Dc = Dc{c, h; e), = Dh{c,h;6) being the general 
form. As far as this paper is concerned, we focus on the simpler dependency cases 
Dc{h; e), Dh{c; e) the extensions being straightforward. For instance, the functional 
forms d?]) become 

D, = l-e^p(^-'^y Dh = I - exp (^-"^y (18) 

where uie^ = Vc/^q, ^2^"^ = I'h/c*^ on taking Rc = h'^, Rh = cP ■ 



3 The sharp interface model derivation 

The limit e — t- will be considered; Ortoleva et al. [27] (cf. sect E, p. 1001) refer 
to this as the fast-aqueous-reaction asymptotics. This corresponds to the bulk re- 
action being very rapid and the reaction is expected therefore to be essentially in 

Natural environments normally contain 0.03—0.05% CO2, with the evolution of the carbonation 
depth being slow and taking many years (typically 5-10 years). An accelerated carbonation chamber, 
exposes concrete to 50% CO2 which dramatically reduces the time needed to perform experiments 
from years to a matter of days (typically upto 20 days). 
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Parameter (units) Value (accelerated) Value (natural) 





= Dl (mVi) 


(0.5 - 5) X 10'^ 


(0.5 - 5) X 10-^ 




(m^s-i) 


10-13 


10-13 




(m s"!) 


1.16 X 10-2 


1.16 X 10-2 


c* 


(mol m-3) 


4.38 


2.71 X 10-3 




(mol m-3) 


1.04 X 103 


1.04 X 103 


k 


(m^mol-^s-i) 


1.74 X 10-5 


4.6 X 10-5 



Table 1: Typical parameter values for natural and acceleratecocarbonation obtained 
from [19], [30], [28], [29]. Concentrations are expressed in units of mol m-3 for 
consistency with the reaction rate (partial orders of the reaction being p = g = 1); 
the molar masses of 44.01 g for CO2 and 74 g for Ca(0H)2 may be used to convert 
to units of g m-3. 



Nondimensional 
parameter 


Value (accelerated) 


Value (natural) 


e2 


10-9 - 10-4 


10-9 - 10-4 




10-11 - 10-^7 


10-1° - 10-5 




10-3 


10-6 


H 


104 - 10^ 


104 - 10"^ 



Table 2: Estimated nondimensional parameters using values in Tabled] 

equilibrium (c = or /i = 0) except in a thin boundary layer, the reaction zone, 
where virtually all the carbonation occurs. In other words, almost all the reaction 
is confined to a very narrow layer, the structure of which will be sensitive to the 
behaviour of the diffusion coefficients. The location of the reaction zone will be de- 
noted by rr = s(t;e). We denote the regions in which the reaction is in equilibrium 
as the outer regions, where outer 1 has /i = and outer 2 has c = 0. The details 
of the reaction zone depend crucially on the behaviour of the diffusion coefficients, 
which we investigate in subsequent sections. For generality in derivation, we keep 
^ = 0(1) and consider its vanishing limit later when seeking to solve the resulting 
sharp interface Stefan problems. 

3.1 Outer solutions 

Taking the limit e — )• and posing regular expansions 

c = co + 0(e2), h = hQ + 0{e^), 
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we obtain two outer regions as below. The location of the reaction zone also requires 
expansion, the leading order term we denote by s{t), the higher order terms of which 
not impacting in the analysis to the orders of the calculations considered. 
Outer 1 < a: < s(t) (carbon dioxide region) 

Here we have ho = (with in fact h=0 to all algebraic orders of e) and 

dco d f dco\ , 

where ~ Dc{co, 0;0). This region is essentially calcium hydroxide free, having 
been used up in the reaction to form calcium carbonate. 
Outer 2 s{t) < x < 1 (calcium hydroxide region) 

In this region we have cq = (and in fact c = to all algebraic orders in e) with 

with D/ig = D/i(0, /iq; 0). This equation may simplify further depending upon the 
relative sizes of 5 and e. This region is almost pure calcium hydroxide, the presence 
of carbon dioxide being negligible. 



3.2 Interface conditions 

In practical terms it is the outer regions that are the ones of most significance; 
however, the interior layer must be analyzed to obtain the continuity conditions 
linking both outers. Without appealing to such an inner analysis, however, we can 
use ([9]) and (fTOj) to obtain the statement 

d , d / ^ dc 6'^ ^ dh\ 

-U.c-f,) = -(^D.---,D,-y (21) 

This may now be used to obtain the jump condition at x = s{t), 

= 0, (22) 
representing conservation of mass. Using the outer solutions this yields, 

atx = s{t), -Dc„ ^ -^^ho = s(/iCo(s ,t) + ho{s^,t)). (23) 

This represents the Stefan type moving boundary condition for ()19p and ()20p . How- 
ever, we still remain two conditions short in terms of specifying the moving boundary 
problem. These final conditions can only be found by undertaking an interior layer 
analysis of the reaction zone. They correspond to kinetic conditions and will take 
the general form 

co{s-,t) = ^i{s), /io(s+,t) = $2(5), (24) 



dc 
dx 
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where the functions and <I>2 are to be determined. To complete the sharp interface 
statement we require suitable initial conditions. These are parameter sensitive, the 
most usual situation being 

at t = 0, Co = Ci{x) for < 2; < s(0), ho = hi{x) for s(0) < x < 1 

and s(0) = sq, (25) 

with, more oftern than not, the initial position of the reaction zone being taken at the 
concrete surface sq = 0- In certain parameter regimes, there is an initial carbonation 
stage in which the reaction zone remains in its initial location for a finite time tQ. 
In the this case we modify to 

at i = Co = Ci{x) for < x < s(0), /iq = Hi{x) for s(0) < x < 1 

and s(0) = sq, (26) 

where Cj and Hi are the resulting concentration profiles at the end of this initial 
carbonation stage. 

The moving boundary problem ([I2D-(12D]), (gSD-dMI) with 1^ or (gU]) is a two- 
phase problem in the regime 5 = 0(e). When 5 <C e, it degenerates to a one-phase 
problem where $2 is no longer needed. 



4 Reaction layer analysis: Slowly varying diffusivities 

We start by considering the simplest case in which neither nor Df^ have any 
explicit dependance upon e. As such they remain order 1, which we refer to as 
slowly varying. The reaction layer in these situations comprises a single region, 
which we describe here for the parameter range (5 < e, larger values of 5 not being 
physically relevant for the carbonation problem. These slowly varying diffusivities 
provide a useful base case against which to compare the effects of more rapidly 
varying diffusivities in later sections [5] and [6j 



4.1 Asymptotic regimes 
4.1.1 Case 5 < eSr 

p+2 

Fullest balance in the governing equations is obtained when 6 = Ae^+i with A = 
0(1). We introduce the scalings 

2 2 
X = s{t) + ep+^x, c = ep+^c, h = h, (27) 

which preserves the CO2 flux from outer 1, to obtain 
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Posing 

c = co + o(l), h = ho + o{l), (30) 

with 

D,{c,h) ^ D,{0,ho), Dh{c,h) Dh(,0,ho), (31) 
we obtain the leading order equations 

The matching conditions for outer 1 and 2 are 

asx^-oo D.^^D,^^^^^^, ho^O, (33) 
ox ox 

dh 

as x — )• +00 Co — )• 0, hQ^hQ(s^,t), X^D^-— — )• 0, (34) 

ox 

and necessarily co(s~,t) = 0. Consequently, we have 

$1 = 0, (35) 

whilst <I>2 is not required. The above analysis also holds for 6 <^ ep+^ where we can 
simply take the limit A — )• 0. 

P+2 

4.1.2 Case ep+i <C 5 <C e 

The scalings in this case are given by 

•e\P_ ^ 

to obtain 



X = s{t) + 6 y-j X, c = ( - ) c, h = h, (36) 



Posing 



"V^/ ^ dx dx \ ^^dx) ' ^'^^^ 

/e\p+i_ - - /e\p+i_ 

c~co + e(^-j ci, /i~/io + e(^^j /ii, (39) 



we obtain the leading order equations 

whilst at first order we have (for brevity we only record the case in which both 
diffusivities constant, additional terms from the expansion of the diffusivities being 
present when they are concentration dependent) 

d [ ^ dci 
dx 

s 



O.^) = P<t"K^i + ifoht'hu (41) 
§^ I (^^> 
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together with the matching conditions for outer 1 and 2 being, 

as X -oo Dc^^O, ^ Dc ^^'^^^ 

ox ox ox 

ho ^0, hi^ 0, (43) 

dh 

as X +00 Co 0, ci —?■ 0, ho ho(s^ ,t), Dh^rz — 

ox 

D,'-^^D,'-h^. (44) 
ox Ox 

Necessarily co(s,t) = and thus ()35p holds with again $2 not required. Here it is 
the first order terms which allow the fluxes to be matched with the outer regions. 

4.1.3 Case 5 = 0(e) 

In this parameter case, an altogether quite different scenario takes place. First there 
is an initial period in which the reaction zone remains at the surface. This period 
may be termed stage I carbonisation, using the terminology introduced for metal 
oxidation [9] where a similar situation occurs. During this stage, the hydroxide at 
the surface is depleted, the details of which are described in Appendix lAl This stage 
ends at a finite time to when the surface hydroxide concentration becomes small 
(i.e. zero at leading order). After this, stage II carbonation takes place where the 
reaction zone moves into the concrete, which we describe here. 
Writing 5 = \e with A = 0(1), the scalings are 



s(t) + gp+'J+ix, c = ep+i+'^c, /i = ep+^+i /i, (45) 



4 2 Qc Q j Qc \ - 

e^'+'+^/i^ -e'^+^+i/is— - — ( Z?c^ j =-0^/19, (46) 

dh 2 _ .dh . 2 / ^ c^/i \ 



givmg 



^^'l^-^^'s—-X'—{Dhj-]=-^h'^. (47) 
ot ox ox \ ox J 

Posing (|3Up . we obtain the leading order equations 

where here Dc{c,h) ~ Dc{0,0), Dh{c,h) ~ 0^(0,0) are constants. The matching 
conditions for outer 1 and 2 being, 

D,^^D,^^^^f^, ho^O, (49) 
ox ox 

- V , - . n n "^^0 . n 9hois+,t) 

asx-^+oo co-^0, Dh^T^ ^ Dc 7: • (50) 

Ox ox 

Now, we have co(s~,t) = = ho{s~^,t) and consequently, 

<I.i = = $2. (51) 
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4.2 Numerical results 



Numerical solutions of the full model equations ([9|)- ()13p are presented here for com- 
parison with the asymptotics. The equations were implemented in the finite ele- 
ment package COMSOL Multiphysics, using the general form PDE mode and a 
mesh of 400 quadratic Lagrange elements (1602 degrees of freedom). The time- 
dependent (Backward Difference Formula) solver had error tolerances abs tol= 10~^, 
rel tollO"'^. For all numerical simulations we take initial data and parameter values 

/i- = 1, Ci = for < x < 1, fj, = 10~^ H = 10^ p = q = 1. (52) 

Figure [2] shows numerical solutions for the case Dc = = 1, in the parame- 
ter regime e = 10"'^. Four selected values of delta have been chosen, covering the 
asymptotic regimes presented in section 14. li There is a clear marked change in be- 
haviour between the case 5 < e shown in (a)-(c) and 6 = e shown in (d). The 
hydroxide profile falls to zero at the reaction front in (d) with noticeable diffusion 
ahead, in contrast to (a)-(c) where it is unity at the front and falls in zero within the 
reaction zone. Moreover, the rates at which the fronts move vary, with a time delay 
apparent in (d) compared to (a)-(c) . This supports the stage I carbonisation regime 
that is postulated to take place when 5 = e, in which the reaction zone remains at 
the surface for finite time as the hydroxide concentration falls. The reaction zone 
starts moving into the concrete once the hydroxide at the surface has been depleted. 




(c) S ^ (d) 5 = e 



Figure 2: Numerical results for slowly varying diffusivities. Tlie parameter values are e = 
10""^, p = q = 1, II = 10^'^ and H = 10'* and selected 6 as stated in (a)-(d). The solid lines 
refer to the concentration of carbon dioxide while the dotted lines are the concentration of 
calcium hydroxide at the shown times. 
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4.3 Sharp interface model summary 

In the case 5 <^ e we obtain the one-phase problem 



dt dx V dx 
dcQ 
dx 

dcQ 
dx 

at t = c = Ci for < x < Si, s = Si, (56) 



mO<x<sit),t>0 ^-± = — {D,,^]; (53) 
onx = -Dc,^ = Hil-co); (54) 

on x = s{t) Co = 0, ~Dcq-t^ = shi{s); (55) 



with Hq = hi for s{t) < x < 1. 

In the case 5 = Xe we obtain the two-phase problem 



inO<x<.(t),t>0 ^-^ = -{D,,^]- (57) 



9co _ / 5co 
dt ~ dx\ '° dx 
dho 2 d ( dho 



i„.(,)<,<M>0 ^ = A^_^«.„^j; (58) 

on I = ""«>^ = ~ 

„ dhr, , , 

onx = l Dh,-^ = 0; (60) 

on X = s(t) Co = 0, ho = 0, "-^co^ = -^^^0^5 (61) 

at t = to Co = Ci for < x < Sj, ho = Hi for Sj < x < 1, 

s = s*. (62) 

Here to is the end of the time at which the reaction zone remains at the outer 
surface and after which it begins to ingress into the concrete. Ci,Hi denote the 
resulting concentration profiles at this time, which differ from their initial values 
Ci , hi respectively. 

5 Reaction layer analysis: A single rapidly varying dif- 
fusivity 

Here we consider the carbon dioxide diffusivity to depend on e as well as the hy- 
droxide concentration, i.e. Dc = Dc{c,h;e). The hydroxide remains slowly varying 
i.e. Dfi = 0(1). The properties required for this diffusivity are 

r 0(e2) iorh = O{l),h>0- 
^ ~ \ 0(1) for h = o(l), h>0. ^^"^^ 

Whilst the functional form is not important, for definiteness we may consider an 
hydroxide only dependent diffusitivity Dc{h;e) in the Arrhenius form 

^c = l-exp('-^], (64) 
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where i^i is a positive constant and Dc{0; e) = 1. [A more general power of h could 
be considered, with the below analysis only needing slight modification with the 
consideration of further terms in the expansions in the inner inner regions.] 



5.1 Asymptotic regimes 
5.1.1 Case <5 < 

We write 6 = Xe^ with A = 0(1). The scalings for the inner inner region are 

X = s{t) + e^x, c = c, h = h, Dc = e^Dc, (65) 

which give 

n dc dc d f - dc\ „-„ , . 

- - ( D,§:) = (67) 

ot OX ox \ ox J 

Posing (|3Up we obtain at leading order 

together with the matching conditions, 

asx-;>-oo cq co(s ,t), D^-^ Dc , no 0, (69) 

dh 

asx^+00 co-^0, ho ho(s~^,t), (70) 

Ox 

An additional region is required to facilitate the matching with outer 1 and 
explain (j69p . due to the order of magnitude change in the carbon dioxide diffusivity. 
We denote this region as inner 1, the scalings for which are 

2 2 ^ 

x = s{t) + eix, c = c, h = eih, (71) 

giving 

We restrict ourselves the parameter case q > 1, the analysis needing slight modifi- 
cation for < g < 1 which is described in Appendix iBl Posing 

2 ^ ^ 

c~co + e9Ci, h ho, (74) 
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we obtain 



^ fr> . .dco d f dcA .dho 



After matching to outer 1 and the inner inner we have 

/ - ,^ ^ r,dci dco{s-,t) - dco 

co = co[s ,t) = $i, = Dc 7. = hm Dc^, 

ox ox x-s>-oo Ox 

: Co{s,tyx 



(76) 



5.1.2 Case < 5 < e 

This case is more involved, where the reaction zone sphts into three regions. In 
progressing from outer 2 to outer 1, we label the regions as inner 2, inner inner and 
finally inner 1 as summarised in Figure [3l 

S ^2^ 



= [CO2] 



h = [Ca(OH)j 



Outer 1 II II a Outer 2 




11 Inner 1 

12 Inner 2 

II Inner Inner 



s{t) ^ X 

Figure 3: Schematic summary of the asymptotic regions for a rapidly varying CO2 difFusivity 
in the parameter case ^ S ^ e. The reaction layer lies between outer regions outer 1 
(0 < X < s{t)) and outer 2 (s(i) < x < 1). It is composed of three regions: an inner inner 

region (II) with 9 = (^^) j together with inners 1 (II) and 2 (12). The inners 1 and 2 
facilitate the matching of the inner inner region with the two outers. The case shown is for 

l + g 

q > 1 and S < Sqr with 5c r = 0(e « ). Minor modifications are needed to inner 1 for the 
cases q < I and S > 6c r- 

The scalings for the inner 2 region are given by 

X = s{t) + (-^^ X, c = 0, h = h, Dc = e^Dc, (77) 
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where A = 0(1) is introduced for convenience and c is actually exponentially small 
in e if p > 1. The governing equation for the hydroxide at leading order in x > is 

s-^ + A — = 0, (78 

ox ox \ ox I 

together with the outer 2 matching conditions, 

dh 

as X — > +00 hQ^ho(s^,t), -rrr- 0. (79) 

ox 



Thus 



au 



and for matching with the inner inner we have the behaviour 

dh 

as X ^ 0+ ho-}- 0, A^Dfe— # s/io(s+, t). (80) 

ox 

In the case Df^ a constant, we have the simple explicit solution 
^0 = ho{s^,t) ( 1 - exp ( --tJt^x 



For the inner inner region, we adopt the scalings 

X = s{t) + (^^^ Bx, c = c, h = eJi, D, = (81) 
where ^ 



— ) « 1- (82) 



The governing equations are now 



At leading order we have 



together with the inner 2 matching conditions, 

asx-^+oo Co ^ 0, Dh^— ^ Dh ^7 , (86) 

ax ox ax 
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and the inner 1 matching behaviour 



asx^-oo co^co(0-,t), ^ 4^^^^, ho ^ 0. (87) 

To ahow matching between the inner inner and outer 1, an inner 1 region is 

required. The scaHngs are determined by the size of 6 relative to a critical value 

(i+g) 

5cr = 0{e 1 ) and the value of q. We consider first the case q > 1, for which 
<^ 6cr ^ £• For 5 < 6cr the scalings ([TT]) apply, as does the subsequent analysis 
culminating in (j76p . This also pertains when 6 = 6cr, with the modification of the 
diffusion term entering the leading order hydroxide equation. For 5 > 6cr, the scalings 
need adjusting to 

f S\ I 1- 
X = s{t) + I — \ ei X, c = c, h = eih, 



Xe 



which give 



Posing 



c ~ Co + ( ^ j e^ci, h^ho, (91) 



we obtain 

' ^A:f) =0, + ^ f l^M) = 0. a4 ( D^-tf] = (92) 



OX \ OX J OX OX \ OX J OX \ ox 

After matching to outer 1 and the inner inner, we again obtain ()76p with, for 
q>l, 

^° - I 2(1 + ,)Z). " J (9') 

in the case when Dh is independent of x. 

For the range g < 1, we note that 6cr < and so it is the regime just discussed 
that occurs with slight modification. Specifically, for q = 1 this inner 1 region is that 
same size as the inner inner and so the spatial scaling in ()88p is modified to 

x = s{t) + -{S{t) + x). (94) 
A 

We again obtain (|75]l ~ ([76]) . but now with 

ho = exp [am , Sit) ^(^^^ log ( (y ) ' 7) , (95) 
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where 

When q < 1, the inner 1 region is narrower than the inner inner. We again have 
()76p and ()93p , where the latter solution is taken for x > and continued by zero for 
X < 0. The inner inner region is now restricted to x > with the inner 1 matching 
conditions (j76]) occurring as x — )• 0+. 



5.1.3 Case 6 = 0(e) 

At the end of stage I carbonation (where the reaction zone remains at the outer 
surface), stage II carbonation takes place in which the reaction zone moves into the 
concrete. 

In this case the reaction zone comprises two regions. We have an inner inner 
region that matches the order 1 flux of hydroxide from the outer 2 region as well as 
allowing the CO2 concentration to fall to o(l). Again an inner 1 region is needed to 
accommodate the fall in CO2 diffusivity, allowing matching between the inner inner 
and outer 1 regions. The details can be deduced from the previous subsection by 

setting 6 = Xe with A = 0(1). Consequently, for the inner inner region the scalings 

2 

(|8ip apply with = e^+i so that again we obtain (j85p as the leading order equations. 
Since there is now no inner 2 region in this case, the matching conditions ()86p now 
apply to outer 2, with ho being replaced with ho. The inner 1 matching conditions 
([87|l remain the same as do the details for the inner 1 region as given by ([88|) - (f95]) . 




Figure 4: A schematic of the asymptotic regions in the case 6 ~ 0{e) of a single rapidly 
varying CO2 diffusivity. The reaction zone is a two layer structure, comprised of an inner 
inner region of width o(^e~^, together with an inner 1 of width O ^e^^ This latter 

transition region accommodates the change in the rapidly varying difFusitivity, the situation 
q > 1 being depicted. 
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5.2 Numerical results 



Here we present numerical simulations for a rapidly varying CO2 diffusivity, using 
the scheme and data/parameter values specified in section IT2l We found it necessary 
to use the following regularisation for the diffusion coefficient ([61 



1 — exp 



+ 91) 



(96) 



with 61 = 10 ^. Although we don't present results for p ^ 1, q ^ 1, we found the 
following regularisation of the reaction term useful to avoid numerical difficulties, 

r = c{c^ + elf-^h{h^ + el)'^, 

particularly for p < 1 or g < 1. 

Figure [5] shows the numerical results for e = lO^'^. Three selected cases of 6 
are shown covering the different asymptotic regimes noted in the previous section. 
Noteworthy features are the 0(1) concentrations of CO2 at the reaction zone in all 
cases. Also the slower progress of the reaction zone compared to the corresponding 
situation in the slowly varying diffusivity case. Again the case 6 = 0(e) shown in 



Figure 5(c) gives different behaviour to the case 5 <C e, with the hydroxide concen- 
tration being small at the reaction zone front and stage I carbonisation taking place 
initially. 




t = 1.0 

t = 2.0 

t = 3.0 

t = 4.0 



0.2 0.4 0.6 0.8 1 



(C) <5: 



Figure 5: Numerical results for one rapidly varying diffusion coefficient Dc as defined in 

for e = 10^^, p = q ~ 1, yU = 10"'^ and H = 10^. The solid line refers to the concentration 

of carbon dioxide while the dotted line is the concentration of calcium hydroxide. 
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5.3 Sharp interface model summary 

In the case 5 < e^, we now obtain the one-phase problem (|53p -(|56 p but with the 
more general moving boundary condition 

on X = s(t) Co = ~^co^ = s(/iCo + hi{s)), (97) 



in place of ()55p . We remark that $i is determined by the solution of the inner inner 
problem 

with 

as X — )• +00 Co — )■ 0, Dc-^ —^0, /lo — ^i(s), (99) 

ox 

that satisfies 

asx^-oo co^^>i(s), /lo^O, 1?^^ ^ 0. (100) 

ox 

For the functional form (|64p we have Z)c = z^i/^o- 

For ^ (5 ^ e we again have a one-phase model with moving boundary condition 
(|97p . but now <I>i is determined by the solution of the inner inner 1 problem 



with 

asx^+cx) co^O, l)c??^0, X^Dh^ ^ shi(s), (102) 

ax ax 

asx^-oo co^$i(i), /lo^O, ^^0. (103) 

ox 

In the case 5 = 0(e), we now obtain the two-phase problem (|57p ~(|62 p but with 
the more general moving boundary condition 

on X = s{t) Co = ^i{s), -Dco^ - ^^^ho-^ = s^co, (104) 

in place of (j6ip . We remark that the outer 2 solution determines $2 = 0, whilst <I>i 
is determined by the solution of the inner inner problem 

with 

asx^+cx) co^O, D,-^^0, ^ By, \ , 106 

ox ox Ox 

dh 

asx->-oo co^^>i(s), 0. (107) 

ox 

This formulation holds for the parameter range q > 1. For q < 1, the boundary 
condition p07p is applied at x = 0"*" and the problem considered on the half-line 
X > 0. 
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6 Reaction layer analysis: Two rapidly varying difFusiv- 
ities 



We now consider the situation in which both the carbon dioxide and hydroxide dif- 
fusivities are rapidly varying. In addition to Dc as defined in section ([5]) we consider 
an hydroxide diffusivity Df^ = Dfi{c; e) with the behaviour 

0. = i2t? !°" = OW'=>°^ (108) 

\ 0(1) for c = o(l),c > 0, ^ ' 

and Dfi(0; e) = 1. This maintains an 0(1) flux of hydroxide in the reaction zone. For 
definiteness, a specific functional form could be the Arrhenius type, 

Dh = I - exp (^-'^y (109) 

where 1^2 is a positive constant. 

6.1 Asymptotic regimes 
6.1.1 Case 6 = 0(e) 

Writing 6 = Xe with A = 0(1), the scalings for the inner inner region are 

x = s{t) + e^x, c = c, h = h, Dc = e^D„ = (^j^ Dh, (110) 

which give 

2 9c .dc d f ^ dc 



dc 


d 


dx 


dx 


i- 


r 


VAe2 





Posing ()30p we obtain at leading order 

.dco d f f, dco\ .dho ^9 f f. dho\ - 

together with the matching conditions, 

D^^^D^^^^^^f^, ho^O, (114) 
dx dx 

as X ^ +00 Co 0, Dh^— Dh ^ , (115) 

dx dx 

Again additional regions are required to facilitate the matching with outer 1 and 
2, due to the order of magnitude change in the diffusivities. As in the single rapidly 
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varying diffusivity case in section ()5.ip we require an inner 1 region, the scalings for 
which are 

2 2. fsV - 

x = s{t)+eix, c = c, h = eih, D^, = I - 1 Dh, (116) 

giving ()72|) -(j73 p , resulting in the same matching conditions (|76|) (with care taken to 
distinguish the cases 0<g<l,g = l and q > 1). 

To facilitate the matching with outer 2, we consider an inner 2 region with the 
scahngs 

x = s(t) + €^{S{t;e)+x), c = epc, h = h, Dc = e^Dc, (117) 
where S(t; e) ^> 1. Thus we have 

We restrict ourselves first to the parameter range p > 1. Posing 

c~co, hr^ho + e^hi, (120) 

we obtain 

^(d^\=0 —(d^\=0 s— + — (d—^=0 (121) 
dx \ ^ dx ] ' dx \ ^ dx j ' dx dx \ dx J 

Thus we have 

f , / 4- X ^ dhi ^ dhQ(s~^,t) - dhn 

ox ox x^+oo Ox 



CO = exp (^--g- j , (122) 

after matching to outer 2 and the inner inner regions, which also determines 

2sb 

S(t;e) ^log(l/e) as e ^ 0. (123) 

P 

For p = 1, slight modification of the above takes place with the leading order reaction 
term entering the equation for cq in (|12ip . For < p < 1, the inner inner solution 
for Co vanishes at x = S{t) say and we modify the spatial scaling in ()117p to 

. = ,w+.^(s«)+.^^ 

As a result the leading order equation for c becomes 

dx V '^^x j ""^^0' 
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with Dc = vi/ho^s^ ,t)'^ which posseses the exphcit solution 



Co = < 



2{i+p)De ^ ) lor X < u, 
for X > 0. 



(124) 



6.1.2 Case 5 <C e 

Only slight modification of the preceeding analysis is required for this case. The 
inner inner region remains the same with the scalings piU|) , but now (|113|) becomes 

dx + dx V ' dx)~ '^''O' ' dx ~ '^''0' ^^^^^ 

with the second condition in (|115p replaced by 

as x — )• +CXD ho^ho{s^,t). (126) 

The inner 1 analysis with scalings (jll6|) remains the same for q > 1, whilst 
for q = 1 the only change is omission of the hydroxide diffusion term which is 
subdominant. There is more of a change for < q < 1 where the spatial scaling 

and balance in the hydroxide equation now depend upon the size of 6 relative to 

(i+g) 

^cr = e 29 . Specifically, 

( 2(1-9) \ 
S{t) + e 9 X ] the dominant balance for the hydroxide 

equation being that in ([75|) for 6 < 6cr and ^ in Appendix [B] when 6 = 6cr', 

(ii) for 5cr < ^ < X = s{t) + ^S{t) + (^)^ e i x^ the dominant balance for 

the hydroxide equation being (fTOj) in Appendix [Bl 

The inner 2 analysis remains unchanged for <C 5 <C e, where X = 6/e <^ in 
()119p . For 5 = e^, the convection term balances the diffusion term in the hydroxide 
equation (|119p . although (|12ip remains unchanged. For 5 <C e^, the hydroxide con- 
vection and diffusion terms balance in a thin region of size x = (5/e^)^ centered at 
S{t), which allows the required matching to outer 2. 



6.2 Numerical results 

We implement the scheme as in the previous numerical sections 14.21 and 15.21 using 
the same initial date and parameter values. An analogous regularisation to (|96p is 
adopted for the hydroxide diffusivity ()109p . 

Figure E] shows the numerical results for e = 10"^ and three selected values of 
5. Compared to the single rapidly varying CO2 diffusivity, the addition of a rapidly 
varying hydroxide diffusivity further slows down the progress of the reaction zone. 
The most significant change occurs in the 6 = e regime, where now the hydroxide 
concentration is 0(1) in the reaction zone and no longer small as for the slowly 
varying and single rapdily varying CO2 cases in sections 14.21 and 15.21 Stage I car- 
bonisation no longer takes place in this case. 
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0.2 0.4 0.6 0.8 1 

X 

(a) 5 = £2 



0.2 0.4 0.6 0.8 1 

X 

(b) 5 = £3/2 




= 2.0 
= 4.0 



0.2 0.4 0.6 0.8 1 



(C) 5 = £ 

Figure 6: Numerical results for two rapidly varying diffusivities. The parameter values are 
£ = 10~2, p = q = 1, fj, = 10^3 and H = 10*. The solid line refers to the concentration of 
carbon dioxide while the dotted line is the concentration of calcium hydroxide. 



6.3 Sharp interface model summary 

In the case 5 = e, we obtain the two-phase problem (|57 p - ()62p but with the more 
general moving boundary condition 

on x = s{t) co = ^>i(s), /lo = $2(s), 

-D^,^-X'D,,^ = s{f^co + ho), (127) 

in place of ()6ip . The solution of the inner inner problem 

■ dco d f - dco\ r,jq .dho d (- dho\ 

with 

as X ^ +00 Co ^ 0, 0, (129) 

ox 

as X — )• —CO ho^O, Dh^^O, (130) 

ox 

determines ^*i,<&2 via 

^>i(s) = lim Co, ^2(5) = lim ho- (131) 

X— 00 X— >+oo 



In the case 5 <C e, we have the one-phase problem (j53j) - (j56|) with (j55|) replaced 

by 

onx = s(t) co = ^i{s), -Z^co^ = s(co + /ii(s)), (132) 
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where <^>i is determined by the solution of the inner inner problem 
with 

as X — +00 Co — )■ 0, ^c-qZ- 0, Tiq ^ hi{s), (134) 
as X — )• —00 Co — )• <I'i(s), TiQ — )• 0. (135) 

For the functional forms ()64p and ()109p we have Dc = 1^1/ Jiq, Dh = 1^2/ 

7 The derived generalised Stefan problems 

7.1 The one-phase problem 

In the case (5 <C e we obtain the one-phase problem 

inO<.<.(t),t>0 ,^ = ^(0,,^); (136) 

onx = -D,,^ = H{l-co); (137) 

onx = s(t) co = $i(s), -L>co^ = s(/iCo + /ii(s)), ; (138) 

at t = c = Cj for < X < Sj, s = Sj, (139) 

with /lo = /li for s{t) < x < 1. For a slowly varying CO2 diffusivity, the analysis of 
section m implies = 0, whilst for a rapidly varying CO2 diffusivity $1 is given as 
the solution of the problem p5|) - (|lU0p or (|lUip - (|lU3p depending upon the size of 5 
relative to e^, or p33p - p35p if the hydroxide diffusivity varies rapidly as well. 

The parameter 1/^ represents the Stefan number, the time scaling t = fir see- 
ing it arise in the Stefan condition in ()138p for its more common occurrence. The 
asymptotics of this problem for power law forms of $1 = are described in [HI 110]. 
Relevant to the concrete situation is the large Stefan number limit — t- 0. In this 
limit, the problem becomes quasi-steady due to the now disparate time-scales of 
diffusion and reaction. This quasi-steady problem holds after an initial transient 
regime t = 0{fi) in which the interface is stationary at leading order. An explicit 
solution for the quasi-steady problem t > 0{fi) when Dc^ = 1 is 

co^^r + ^l^, (140) 

with the interface s(t) determined by the solution of 

hi{s)(^ + s]s + <^>i{s) = l, s{0) = s^. (141) 
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7.2 The two-phase problem 

In the case 5 = Ae we obtain the two-phase problem 



dco d / dc{ 



inO<.<»«).<>0 „rni = __^^z,^,r£j; (142) 
in.W<.<M>0 (143) 

onrE = -D^^ = H{l-coy, (144) 

onx = l Dho^ = 0; (145) 

onx = s(t) co = ^>i(s), /io = $2(s), 

-D^,^-X'D,,^ = s{f^co + ho); (146) 
at t = to Co = Cj for < x < Sj, /iq = i/j for Sj < x < 1, 

s = Si. (147) 

Here to is the end of the time at which the reaction zone remains at the outer 
surface and after which it begins to ingress into the concrete. Ci,Hi denote the 
resulting concentration profiles at this time, which differ from their initial values 
Ci,hi respectively. For slowly varying CO2 and Ca(0H)2 diffusivities, the analysis 
of section m implies = = $2- For a single rapidly varying CO2 diffusivity 
is given as the solution of the problem ()105p - ()107p with $2 = 0. In both these 
situations there is a stage I carbonisation period where the Ca(0H)2 concentration 
in the initial reaction zone at the surface falls to zero. When both diffusivities are 
rapidly varying then $1, $2 are determined by ()128p - ()130p . there not being a stage 
I carbonisation period in this case. 

In the quasi-steady CO2 limit /U — )• 0, we again have solution ()140p in the case 
Dco = ^- This gives the non-standard Stefan problem for the Ca(0H)2 as 

m.(^)<x<l,t>0 ^ = A^^(^.o^J; (148) 

onx = l Dh,^ = 0; (149) 




at t = to ho = Hi for Sj < X < 1, s = Sj. (151) 



This problem in the case slowly varying diffusivity case when = = ^2 has 
been considered by [11] in the context of binary alloy oxidation. However, the more 
general statement has not received attention. 

8 Discussion 

A set of reaction-diffusion equations representing the concrete carbonation process 
has been analysed in the limit of fast reaction - slow diffusion. Different sharp in- 
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terface models in the form of generalised Stefan problems are derived, depending 
upon the properties of the diffusivities (of the two main species CO2 and Ca(0H)2) 
as well as the size of the relative transport parameter = D^hP /D^C* < 0{1). 

These two considerations determine the type of sharp interface Stefan and kinetic 
conditions. A one-phase model results for 5 <C e and a two-phase model when 
5 = 0(e). The sharp interface kinetic conditions are determined by the reaction- 
diffusion equations within the thin reaction zone, thus capturing features of the 
model on smaller length scales (microscales) . The resulting sharp interface model 
considered as a macroscale model thus contains microscale information, which we 
argue is more appropriate than simply stating an empirical kinetic condition. 

The issue of which diffusivity situation is appropriate i.e. one or both rapidly 
varying, is a modelling aspect. The choice will ultimately depend on the best fit 
with experimental observations for the carbonation zone. 

It is worth noting the following open issues: 

• In |15j (sect. 1.5 and chapter 9), one derives a class of distributed-microstructure 
models. Roughly speaking, it is about two-scale systems of PDEs which are 
coupled via micro-macro boundary conditions. [T7] and [18] report on such 
two-scale problems where, at the micro level, a sharp-interface model is re- 
sponsible for the evolution of micro-free reaction boundaries. Can one derive 
via homogenization techniques our two-scale model? 

• Are there any connections between our two-scale result and the Localized 
Model-Upscaling (LMU) method by Pierre Degond et al. [7]? LMU is typically 
used in the context of Boltzmann-like equations. Essentially, it consists in cou- 
pling a perturbation model and its asymptotic limit model via a transition 
zone. In the transition zone, the solution is decomposed into a microscopic 
fraction (described by the perturbation problem) and a macroscopic one (de- 
scribed by the limit problem). 

• The single scale Stefan- like problems can be investigated by means of standard 
methods for free-boundary problems; see e.g. [11]. However, the well-posedness 
of the two-scale FBPs summarised in Section [7] is an open problem. We expect 
that methods developed in [21 [2ll [22] and [3] can be adapted to cope with the 
two-scale structure of our problem. 

• Once the limit problems are shown to be well-posed, the next step is to prove 
rigorously the passing to the limit e — by getting an upper bound on the 
convergence rate (corrector estimate) . This issue is open even for the derivation 
(as fast-reaction limit) of the standard Stefan problem. 

• For the practical corrosion problem, the geometry of the concrete structure 
plays an important role. We would like to extend the asymptotics to 2D and 
3D cases and explore the effect of corners on the speed of the moving reaction 
interface. 

• Resolving numerically in an accurate manner 2D and 3D singular-perturbation 
scenarios require a very good understanding of the evolution of the singular 
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part - the reaction layer. We hope that the asymptotic method will help us 
to construct a good method to capture the a priori unknown position of the 
layer. 

• Using eventually the same experimental data as reported in or [25], we 
want to validate the two-scale model. It is unclear for the moment what is 
the best option: Should one use the two-scale sharp interface model (which is 
exact), or its e-approximation (which is right only within a corrector range to 
be established)? 
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Appendix A. Stage I carbonisation 
A.l Slowly varying diffusivities 

Here we complete the analysis in the case 5 = 0(e) considered in section 14.1.31 
for slowly varying diffusivities. The reaction zone remains at the surface until the 
hydroxide within it is consumed. As such h = 0(1) and for this inner region we 
consider the scalings 



2 



2 



X = X 



c = ep+i c, 



h = h 



(1) 



givmg 




(3) 



Posing 



2 



h = ho + ep+i/ii H 



2 



c = Co + ep+1 ci + • • • 



(4) 



31 



after an initial transient t = 0(ep+i ), we obtain at the leading order 



dx V dx 



dx 



dx 



subject to 



at X = 
as — 7- +00 



ox ox 
Co 0, ho^ho{0,t). 
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(7) 



Consequently, 



ho = ho{0,t), Co 



^ (1 + ^a,it)x 
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D. 



^ exp {-ai{t)x) , 
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(l-p) > 



(i-p) 



p > 1, 
p = l, 

< J9 < 1, 
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where 



2/io(0,t)'? 



[l+p)D, \D, 



p— i\ i+p 



and setting p = 1 gives cri(t). This solution recorded for cq being relevant when Dc 
is constant in this region. 

At the next order for the hydroxide, we have 



with 



giving 



at X = 



ox 



dx 



dhi 

= 0, 

ox 
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dx 
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The leading order outer 2 problem is thus the fixed domain problem 
inO<x<l,t>0 

Ot O2 
atx = x^Dh^ = H; 

OX 



— = X^— ( Dh— 

dt dx \ dx 



at X = 1 
at t = 



dx 

/lo = /ii for < X < 1, 
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the condition ()13p following from matching with the above two term inner solution 
dl]) and using pTj) . In the case when and hi are constants, we have the explicit 
solution 



HI \ 

ho{x, t) = hi-Ht + [^{2-x) + Y,fn eM->''^Dhn\h) cos{mTx) , (16) 

n=0 / 



h 



where the fourier coefficients are 



1 2 

/o = i:, fn = -^r^ for n>l. (17) 

In the semi-infinite case, the corresponding solution is 

h ( ,\ h ^ /■* exp(-xV4A^A^^ 
ho{x,t) = hi Y r= 

We denote by t = to the time at which the leading order hydroxide concentration 
falls to zero at the reaction zone i.e. /io(0,t) = 0, the resulting concentration profile 
then being Hi{x) = ho{x,tQ). In the finite domain case, p6|) gives a transcendental 
equation for to, 

H °° 

HtQ = K + -J— V /„ exp{-X^Dhn\\), (19) 

n=0 

with then 

Hi{x) = '^i-Hto+^^ i-x{2 - x) + Y,Jnew{->''^Dhn^Tr\)cos{mTx)j . (20) 
In the semi-infinite domain case, ()18p gives explicitly 

*o = MJ2 ' -^^a;) = r / ^= ds. 21) 

Since the leading order inner solution for the hydroxide Iiq does not vary spatially 
within the inner region, this Stage I carbonation persists until the leading order outer 
2 solution falls to zero at the reaction zone. Subsequently Stage II carbonation takes 
place as described in section H. 1.31 

The case H = oo \s different from the above, since c as well as h are now order 1 
in the reaction zone. 



A. 2 Single rapidly varying diffusivity D^. 

For the single rapidly varying CO2 diffusivity considered in section I5.1.3t a stage 
I carbonation regime takes place analogous to the slowly varying case. The CO2 
diffusivity is now small in the reaction zone, consequently the scalings ([1]) change to 

4 4 _ _ 

X = e^x, c = e^c, h = h, Dc = e^Dc, (22) 
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and the expansion modified to 

4 _ _ 4 _ 

c = co + ep+iciH , /i = /lo + ef+i/ii H . (23) 

The rest of the analysis remains similar, with Dc replaced by Dc which we note 
for the specific form (|64p is = i'i/h(){0~^ ,ty . This pertains until /io(0+,t) falls 
0(e^/'^), the waiting time to being the same in this case. 



Appendix B. The inner 1 asymptotics 

We complete here the analysis of the inner 1 region in section 15.1.11 for the pa- 
rameter regime < g < 1. 

In the case q = 1, this inner 1 region is the same size as the inner inner region, 
although located in its far-field. We scale using 

X = s{t) + e^{S(t) + x), c = c, h = e^h, (1) 

where 1 <C S{t) <^ 1/e^. The governing equations become 

Posing ()74p with g' = 1, we obtain ()75p with the equation for /iq being 

After matching to outer 1 and the inner inner we again obtain (|76p with /iq and S{t) 
(in the case Dh constant) given by 

2 

ho = exp {a{t)x) , S{t) r — log(e), (5) 

qa[t) 



where 



-s + {s^ + 4X^DhCo{s~,ty)"2 
a{t) = 



In the case < g < 1, the inner 1 region is smaller than the inner inner. We scale 
using 

X = s{t) + {S{t) + e~T- , c = c, h = e^h, (6) 

where S{t) < can be taken as the location at which the inner inner solution Iiq 
first becomes zero. We thus obtain 

idh i^,. ..dh o 9 / dh\ 
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Posing 

(i+g) . . 

c~co + e 9 ci, hr^ho, (9) 



we again obtain ()75p with the equation for /iq now being 

Matching to outer 1 and the inner inner gives ()76p with /iq (in the case Dh constant) 
given by 

ho={ CV(f;glr.f ^^)^ forf>0, (11) 
D for a; < 0. 
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